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We consider the scenario where a 4- lattice constant, rotationally symmetric charge density wave 
(CDW) is present in the underdoped cuprates. We prove a theorem that puts strong constraint on 
the possible form factor of such a CDW. We demonstrate, within mean-field theory, that a particular 
form factor within the allowed class describes the angle-resolved photoemission and scan tunneling 
spectroscopy well. We conjecture that the "large pseudogap" in cuprates is the consequence of this 
type of charge density wave. 



I. INTRODUCTION 

After almost two decades of experimental study, it 
is known that the high temperature superconductors 
have the foUowing known ordered states: 1) antifer- 
romagnetic order at very low doping (x ^ 3%), 2) 
the d-wave superconducting (DSC) order for 5% < 
X ^ 30%. While these two orders exist in all 
families of cuprates, there is a third order, namely, 
3) a 4-lattice constant charge and 8-lattice constant 
spin density wave order, occurring near doping x = 
1/8 in the Lai.48Ndo.4Sro.i2Cu04/Lai.875Bao.i25Cu04 
(LNSCO/LBCO) systems^ There is a wide-spread belief 
that this charge/spin density wave order is anisotropic, 
i.e., they form stripes^^^. 

A significant part of the high-Tc mystery lies in the 
behavior of the underdoped systems'^. Based on specific 
heat^, nuclear magnetic resonance^ , DC transport^, op- 
tical and Raman spectroscopy^' angle- resolved photoe- 
mission (ARPES)^^ and tunnehng^^"^^ , Tallon and Lo- 
ram have made the case that the high-Tj, superconductors 
possess two energy gaps, a pseudogap and a supercon- 
ducting gap^. Recently ARPES experiments on LSCO 
systems^^ and underdoped Bi2212^^ both point to a large 
pseudogap in the antinodal region and a superconduct- 
ing gap near the Brillouin zone diagonals. Similar result 
has also been found in electronic Raman scattering ex- 
periment on Hgl201^^. In addition, it is shown that for 
underdoped Bi2212 a large pseudogap exists in the antin- 
odal region even at temperature w 3Tc, while a gapless 
Fermi arc exists near the nodes. 

Recently, there are clear evidences from the the scan 
tunnehng spectroscopy (STM) studies suggesting the 
presence of a 4-lattice constant checkerboard order in 
NaCCOCi" and underdoped Bi2212i5'i6. Interestingly 
ARPES study has shown that in Na2:Ca2-a;Cu02Cl2, 
where STM found checkerboard order the Fermi arcs 
survive^^. 

In view of these new experimental results we ask the 
question "can the pseudogap in underdoped cuprates be 
caused by some kind of checkerboard CDW?" . To answer 
the question, we will look at the effects of the checker- 



board CDW on low energy quasiparticles. Since the ex- 
istence of low energy quasiparticles is an experimental 
fact, it is reasonable to model the influence of CDW by 
an effective scattering Hamiltonian of the form 

HcDW = ^5]^[/(k,Q)C+,Q^Ck. + /i.c.], (1) 

Q k cr 

where Q is the CDW ordering wavevector and /(k, Q) 
the form factor. 

In the following, we will flrst explore the symmetry 
property of the checkerboard CDW form factor using 
the experimentally observed STM patterns in Sec. II. 
In Sec. Ill, we compare the low energy ARPES and 
STM spectral functions generated by two representatives 
among the allowed form factors. Section IV is the sum- 
mary. 

II. TWO THEOREMS ABOUT /(k, Q) 

In Fig. 1(a), we reproduce the STM dl/dV image of 
NaxCa2-a;Cu02Cl2 from Ref.l4. This particular image is 
made at bias voltage 30 mV. However, the same checker- 
board pattern was seen in a wide bias range — 150mV < 
V < 150mV. Experimentally, it was determined that 
such a checkerboard pattern contains ±Q, where 

Q = (2V4,0),(0,27r/4) (2) 

as its fundamental ordering wavevector. Hence we limit 
the Q summation in Eq. (1) to those given by Eq. (2) 
and k to the first Brillouin zone. In Fig. 1(b), we repro- 
duce the two-point correlation function of the observed 
image presented in Ref.l4. From Fig.l(a,b) we construct 
a caricature in Fig. 1(c) to capture the essence of the 
observed checkerboard. Interestingly, in each 4x4 unit 
cell there are two inequivalent centers about which the 
checkerboard is symmetric under 

Civ = {E, 6*2, 

}, (3) 

the point group of the square lattice. (Here E represents 
identity, and C2,4 denote 180 and 90 degree rotations. 
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FIG. 1: (Color online) STM dl/dV map (a) and the auto- 
correlation image of < lOOmeV LDOS maps (b) from 
Hanaguri et al}^ on NaCCOC, showing the 4x4 ordering, 
(c) Caricature of the observed image shown in (a), (d) Possi- 
ble LDOS pattern which exhibits 6 independent intensities in 
the 4x4 unit cell. In panels (c) and (d), two nonequivalent 
s-symmetry centers are indicated by arrows. In panel (d), the 
d-symmetry centers are indicated by the ellipses. 



and a denotes reflection.) In the following, we take this 
as implying that Hqdw is C4^-invariant about these two 
centers. 

Theorem I A CDW that has the ordering wavevectors 
given by Eq. (2) and possesses a center of Civ symmetry 
in its unit cell must have the following properties. 

(1) There must exist another inequivalent C^y center 
in the unit cell. This second center is displaced from the 
first by the (2,2) translation or its equivalent. About 
these two centers /(k, Q) has s-symmetry. 

(2) There must exist two other centers around which 
HcDW remains invariant under C2v, the subgroup 
formed by the first four elements of Ci^, but changes 
sign under CijC^jax+y^cTx-y Spatially these two 
new centers must be displaced from the two C^^ 
centers by the (2,0) and (0,2) translation or their equiv- 
alents. About these two centers /(k, Q) has d-symmetry. 

Proof. Let us assume Hqbw is invariant under C'^y at 
the origin, i.e., TZHccmVr'^ = Hqdw where TZ G Ciy. 
This implies 

/(7^k,7^Q) = /(k,Q) (4) 

from Eq.(l). After a translation t the form factor changes 
to 

/(k,Q)^5(k,Q) = /(k,Q)e^Q-*, (5) 



where t can be any one of the 16 possible displacements 
within the unit cell 

t — {171,71), TO, n = 0,1,2,3. (6) 

Due to the fact that Q only takes one of the four possible 
values given in Eq. (2) it can be easily checked that for 
t = (2,2) 

ff(7^k,7^Q) = g(k,Q) ynec^y, (7) 

and for t = (2,0), (0,2) 

gink, TZQ) .g(k, Q) K e C2v 

g(7^k,7^Q) = -g(k,Q) Tie 0^,-02., (8) 

where C^^ - C2v = {C4, Cf, (Jx+y, (Jx-y}- QED 

Theorem I implies that any 90-degree rotationally 
symmetric CDW with (±27r/4, 0), (0, ±27r/4) ordering 
wavevectors must simultaneously possess s-symmetry 
centers and d-symmetry centers. The presence of both 
symmetry centers is a necessary consequence of the 
CDW being rotationally symmetric. This fact was 
overlooked in the earlier version of this paper. Con- 
versely any four lattice constant CDW that does not 
possess both symmetry centers must break rotation 
symmetry. In addition, it can be shown easily that a 
rotationally symmetric CDW discussed above possesses 
6 inequivalent sites in the unit celP^, hence allowing 6 
different values of dl/dV. This is shown in Fig. 1(d). 

Theorem II If i/cDW is time reversal invariant, 
/(k, Q) must be real if one chooses either d- or s- 
symmetry center as the origin. 

Proof. Time reversal symmetry requires 

r(k,Q)-/(-k,-Q). (9) 

Since /(k, Q) is invariant under the 180 degree rotation 
about the s and d centers we have 

/(k,Q) =/(-k,-Q). (10) 

As a result, 

r(k,Q) = /(k,Q), (11) 

i.e, /(k, Q) is real. QED 

III. EFFECTS OF THE CDW ON ARPES AND 
STM SPECTRA 

In this section we apply the two theorems proven above 
and take the input from a previous renormalization group 
calculation^"^ to guess the plausible form of /(k, Q). We 
then investigate the effect of the checkerboard CDW on 
the STM and ARPES spectral functions of the low en- 
ergy quasiparticles. We stress that the purpose of this 
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section is not to prove that the ground state of certain 
microscopic Hamiltonian has CDW order. Rather, we 
take a phenomenological approach by assuming its exis- 
tence and look at its consequences that are observable by 
STM and ARPES. 

In Ref.23 it was shown that, with the help of electron- 
phonon interaction, a class of electron-electron scatter- 
ing is enhanced at low energies. This class of scatter- 
ing involves (momentum conserving) scattering of a pair 
of quasiparticles near the antinodes. For example, con- 
sider a pair of quasiparticles lying on the opposite sides 
of the almost nested Fermi surface near the (tt, 0) antin- 
odes as shown in Fig. 2(a). After the scattering these two 
quasiparticles switch sides. The momentum transfer in 
such a scattering is the "nesting wavevector" of the antin- 
odes. For systems such as NaCCOC^^ and underdoped 
Bi2212^^ it has been shown that such nesting wavevec- 
tors are approximately given by Eq. (2). Interestingly, 
Ref.23 also shows that accompanying each such scatter- 
ing there is a related process, whose scattering amplitude 
has opposite sign, where one of the quasiparticle scat- 
tering takes place near the (0, tt) rather than the (tt, 0) 
antinode [Fig. 2(b)]. It was also noticed that when this 
type of quasiparticle scattering grows strong it tends to 
drive a CDW whose form factor has the property that 

sign[/(7^k, Q)] = -sign[/(k, Q)] for 7^ e - ^2,(12) 

In the following, let us choose the rf-symmetry center 
as the origin. Thus 

/(7^k, HQ) = -/(k, Q) for UeC^^- C2v. (13) 

Combine Eq. (13) with Eq. (12) we obtain 

sign[/(k, 7^Q)] = sign[/(k, Q)] for 7^ e - (14) 

In addition, Eq. (12) plus the continuity condition re- 
quires 



/(k, Q) = for k along x±y. 



(15) 



The above considerations lead us to the following ansaz 
for the CDW form factor 

/(k, Q) = 5k(Q)(cos fc^ - cosky) = 5k(Q)/o(k), (16) 

where 6'k(Q) > 0. In the following we shall pick a simple 
realization of Eq. (16) and focus on k lying close to the 
Fermi surface. 

In general the CDW couples each k to other 15 k points 
in the first Brillouin zone. However, most of these 16 
k's lie far away from the Fermi surface, hence can be 
omitted in the low-energy theory. This suggests that one 
only needs to keep a few close neighbors for each k. An- 
other important consideration guiding our construction 
of HcBW is the requirement that a robust antinodal CDW 
gap exists for reasonable change of doping. It turns out 
that this requirement is satisfied as long as the nested 
scattering across the antinodal Fermi surface is the dom- 
inant scattering process. 
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FIG. 2: The two enhanced sets of electron-electron scattering 
(panels (a) and (b)), as obtained from a renormalization group 
calculation^^. The scattering amplitude between these two 
sets differs by a sign, (c) The CDW-inducod quasiparticle 
scatterings (only those in the first quadrant of the Brillouin 
zone arc shown) . The solid lines in these figures represent the 
normal state Fermi surface. 



Put all the constraints together we consider the follow- 
ing quasiparticle Hamiltonian in the absence of supercon- 
ducting pairing 



W = ^*+(k)^(k)*<.(k), 



(17) 



k,cr 



where 



(k) = c++Q^_,, c+_Q^^ J, (18) 



and 



^(k) 

In Eq. (19) 



/ ek 5o/o(k) 5i/o(k) 52/o(k)- 

<5o./'o(k) Ek+Q, 

5i/o(k) ek+Q. 

\52/o(k) ek-Q, 



.(19) 



Qi = -sign(fc,)(2V4,0),Q2 = (0,2^4) for |fc,| < \ky\ 
Qi = -sign(fcy)(0,27r/4),Q2 = (27r/4,0) for |fc,| > \ky\ 

as shown schematically in Fig. 2(c). In addi- 
tion, we expect So to be stronger than Si and 

^2. For the normal state dispersion, we use 
Ck = + ti[cos{kx) + cos(fcy)]/2 -I- ^2 cos(fca;) cos(fcy) -I- 
t3[cos{2kx) + cos(2fcy)]/2 -I- t4[cos{2kx) cos{ky) + 
cos{2ky) cos{kx)]/2 + t^ cos{2kx) cos{2ky) , with 
the hopping constants (in eV) {ti,...,t-^) = 
(-0.5951,0.1636,-0.0519,-0.1117,0.0510)25. In 
the following, we will compare the effects of the CDW 
for the two cases where the Fermi surface is nested/not 
nested by the Q given by Eq. (2). (We adjust to to 
control the degree of nesting.) As to the CDW order 
parameter, we choose 



So = Ac, Si = sAc, S2 = sAc. 



(20) 



We first discuss the case with Fermi surface nesting. In 
Fig. 3 (a) and (c) we present the real space dl/dV image at 
bias voltage 20 mV and the ARPES intensity map at the 
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FIG. 3: Panel (a) and (b) arc the dl/dV images for the CDW state where the /()(k) in Eq. (17) and Eq. (19) is cos kx — cos ky 
and \coskx — cosfcy|, respectively. The window of view is 15 x 15 lattice unit cells. Panel (c) is their Fermi energy ARPES 
intensity maps (Both form factors give the same intensity map). Here only the first quadrant of the Brillouin zone is shown. 
In making these figures we have chosen Ac to produce a 60 meV gap at the antinodes. The parameter s in Eq. (20) is chosen 
to be 0.2. A quasiparticle energy broadening of 10 meV and a to = 0.0945 eV are used. 



Fermi level. These results are calculated with s = 0.2 in 
Eq. (20). The primary effect of changing s is to 1) change 
the intensity variation in the black perimeter in each unit 
cell in Fig. 3(a); and 2) affect the strength of shadow 
band in Fig. 3(c) (see later). Except these changes, the 
main features of both results are preserved. In Fig. 3(b) 
we show the dl/dV image resulting from Eq. (17) where 
the /o(k) in Eq. (19) is replaced by | cos /c^ — cos/cyj (Of 
course, after such a choice the d-symmetry center be- 
comes the s-symmetry center). The purpose of this figure 
is to demonstrate the sensitivity of the real space image 
on the sign of /q. Indeed, while the ARPES image is 
completely unaffected by such a change, the real space 
dl/dV is strongly modified. Upon a comparison with the 
checkerboard pattern observed in Naa;Ca2-£!;Cu02Cl2^'*, 
it is clear that the form factor cosfcx — cos ky (Fig. 3 (a)) 
produces the real space description best. To better un- 
derstand the Fermi arc present in Fig. 3(c) we note that 
in the presence of CDW, the new Fermi surface is deter- 
mined by 



det[A{-k)] = 0. 



(21) 



Since fie<[^(k)] is real (because A(k) is Hermitian) 
dei[^(k)] = yields a single equation with two unkowns 
{kx and ky). Gencrically, one expects the solutions to 
form c/osed one-dimensional curves. Since /o(k) vanishes 
at the node, it is natural to expect the Fermi surface to 
be practically unaffected in its vicinity. Such an unaf- 
fected piece of the Fermi surface and its CDW shadows 
form a closed contour. The reason that in Fig. 3(c) only a 
Fermi arc is visible is due to the CDW coherence factor^^. 
In Fig. 3(c), the strongest shadow band effect shows up 
near the end of the Fermi arcs. Note that such shadow 
band position is very different from that expected from 
antiferromagnetism. Presently there is no report of see- 
ing such shadow bands^"'^^. The reason may be: 1) the 
CDW correlation length as observed by STM experiment 
is not sufficiently long (It is typically of 10 nanometers); 
2) in the pseudogap regime, the superconducting pairing 
still persists. In all cases we studied, the superconducting 



pairing is very effective in weakening the shadow band ef- 
fect. When moving away from the zero binding energy, 
we find that the main changes in the ARPES intensity 
map are: 1) the intensity in the antinodal regions in- 
creases, and 2) the Fermi arcs shrink and move towards 
the origin of the first Brillouin zone. 

By considering all panels of Fig. 3, it is obvious that it is 
the checkerboard CDW with /o(k) = cos kx — cos ky that 
reproduces both the ARPES and STM phenomenology 
well. Therefore, we will only consider this kind of form 
factor in the rest of the paper. 

Now, we turn to the case without Fermi surface nest- 
ing. In this case, using the checkerboard CDW with an 
order parameter of the same magnitude as that in Fig. 3, 
we obtain a weaker fragmentation of the Fermi surface 
as shown in Fig. 4(a). As to the real space pattern (not 
shown), the only difference with Fig. 3(b) is a slight in- 
crease in the intensity variation in the dark perimeter 
region. 

Next, we turn on a DSC pairing and ask what is the 

signature of the checkerboard CDW and superconducting 
pairing coexistence in STM. In this case the Hamiltonian 
becomes 



where 



and 



H = ^<i>+(k)H(k)a>(k), 



$+(k) = («'+(k),vl/^(-k)), 



_ / ^(k) B{k) 
- I B*{\<i) -A(-k) ; • 



(22) 



(23) 



(24) 



In the above equations 
By(k)=Ofori^j 

-Bii(k) = Ak, Ak+Qi, Ak+Qs, Ak-Q2 for i = 1,2,3,4. 

(25) 
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FIG. 4: (a) The ARPES intensity at Ef for the checkerboard CDW with the form factor cos k:r — cos ky and a 60 nieV gap. (b) 
and (c) are the q = (0, 0) and q = {n/2, 0) Fourier components of LDOS for a state with a 60 meV CDW gap and a Ao = 40 
meV DSC pairing parameter. A to = 0.1215 eV is used. 



For d-wave superconducting (DSC) pairing Ak = 
Ao(cos kx — cos ky)/2. In the presence of inversion sym- 
metry (^(— k) = A(k)) the Hamiltonian in Eq. (24) can 
also be written as 

W(k) = A(k) (g) CT3 + B(k) (g) CTi . (26) 

In that case because W(k) anticommutes with I ®ay the 
cigcn spectrum is particlc-holc symmc;tric,. Under stich 
condition the zero-energy eigenvectors are also eigenvec- 
tors oi I ® ay. As a result, the locus of zero energy 
satisfies 

det[A{W, ± iB(k)] = 0. (27) 

Since this determinant is complex, setting its real and 
imaginary parts to zero gives two equations for the two 
unknown kx and ky. Consequently, one expects the so- 
lutions to be isolated points in the Brillouin zone. Thus 
with the DSC pairing the Fermi arc produced by checker- 
board CDW is reduced to point gap nodes. 

In Fig. 4(b) and (c) we consider the case where a 60 
meV checkerboard CDW order parameter coexists with 
a Ao = 40 meV DSC pairing. Fig. 4(b) shows the spatial 
averaged local density of states (LDOS). Note that the 
CDW feature on the negative bias side is much weaker 
than that of the positive side. This is because it is over- 
whelmed by the density of states due to the van Hove sin- 
gularity. The two peaks on the positive bias side are the 
original antinodal coherence peak split by the CDW or- 
der. We have checked that the energy separation between 
these peaks is proportional to the CDW order parameter. 
Another way to determine the strength of the CDW or- 
der is to Fourier transform LDOS at the CDW ordering 
wavevector. In Fig. 4(c), the real part of the q = (7r/2, 0) 
component of LDOS is shown. The two peaks on the pos- 
itive bias side of Fig. 4(b) now appear as a peak and an 
anti-peak. Again, the distance between them is propor- 
tional to the CDW order parameter. Thus we propose 
that by studying the Fourier transformed LDOS, it is 
possible to extract the strength of CDW ordering. 

In Fig. 5 (a), we show several ARPES momentum 
distribution curves (MDC) along the momentum cut 




FIG. 5: (a) The ARPES MDC along the momentum cut 
(— 7r/2,7r) — »■ (7r/2,7r) at three different energies in the pres- 
ence of a checkerboard CDW state, (b) The energy gap along 
the normal state Fermi surface for the pure checkerboard 

CDW state (dashed curve) and a state with coexisting CDW 
and DSC order (solid curve). The form factor of the CDW 
is cos kx — cos ky and the CDW parameters are the same as 
those in Fig.3. 



(-7r/2,7r) ^ (7r/2,7r) for the checkerboard CDW. All 
energies considered here are below the CDW gap. The 
presence of two non-dispersive MDC peaks separated by 
the CDW ordering wavevector is apparent. This is very 
similar to that observed in Ref.21. 

In Fig. 5(b), we present the energy gap along the nor- 
mal state Fermi surface for a pure checkerboard CDW 
state (dashed curve) and a state with both checkerboard 
CDW and DSC pairing (solid curve). The purpose of 
this figure is to illustrate the effect of DSC pairing in the 
pseudogap state. It shows how Fermi arc is replaced by a 
gap node. With thermal phase fluctuations, this explains 
why Fermi arcs shrink to four points as temperature ap- 
proaches zero as observed recently^*^. Given these results, 
we feel quite tempted to associate the larger checkerboard 
CDW gap with the large pseudogap and the smaller pair- 
ing gap on the Fermi arc with the small pseudogap. 

In the literature it is widely believed that the pseudo- 
gap is a consequence of the short-range antiferromagnetic 
correlation'\ Thus it is natural to ask what is the rela- 
tion between the checkerboard CDW discussed above and 
such physics. On microscopic level the CDW presented in 
this paper represents the modulation in the hopping (or 
antiferromagnetic exchange) integrals. Consequently, it 
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is a kind of spin Pcirls distortion which, of course, is com- 
patible with the spin singlet pairing tendency of a quan- 
tum antiferromagnet. In addition to the above remarks 
we note that in a recent paper^^ it is found that checker- 
board CDW is a self-consistent solution of a t — J like 
model at mean-field level, again testify that checkerboard 
CDW does not contradict the superexchange physics. 

IV. CONCLUSION 

In this paper, we present a symmetry constraint on 
the form factor of a 90 degree rotationally symmetric, 
commensurate, checkerboard charge density wave. Fur- 
ther guided by a previous renormalization group study^^ 
we construct a simple model describing the scattering 
of the low energy quasiparticlcs by the CDW. We then 
calculate the low energy ARPES and STM spectra us- 
ing this simple model. The results compare favorably 
with the existing experiments. In particular, the re- 
sults show a spatial dl/dV pattern similar to the one ob- 



served in Naa;Ca2-2:Cu02Cl2 and underdoped Bi2212 by 
grpjyj-i4,i6_ Moreover, in the momentum space it produces 
Fermi arcs resembling those observed by ARPES^^'^^. In 
the presence of a d-wave superconducting pairing, the 
Fermi arcs of the checkerboard CDW are reduced to 
four gap nodes^". Therefore, this study supports the no- 
tion that the large antinodal pseudogap in underdoped 
cuprates is generated by the checkerboard charge density 
wave^^^^ conjectured at the beginning of the paper. 
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